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We consider the passage of long polymers of length N through a hole in a membrane. If the process 
is slow, it is in principle possible to focus on the dynamics of the number of monomers s on one side 
of the membrane, assuming that the two segments are in equilibrium. The dynamics of s(t) in such 
a limit would be diffusive, with a mean translocation time scaling as N 2 in the absence of a force, 
and proportional to N when a force is applied. We demonstrate that the assumption of equilibrium 
must break down for sufficiently long polymers (more easily when forced) , and provide lower bounds 
for the translocation time by comparison to unimpeded motion of the polymer. These lower bounds 
exceed the time scales calculated on the basis of equilibrium, and point to anomalous (sub-diffusive) 
character of translocation dynamics. This is explicitly verified by numerical simulations of the 
unforced translocation of a self-avoiding polymer. Forced translocation times are shown to strongly 
depend on the method by which the force is applied. In particular, pulling the polymer by the 
end leads to much longer times than when a chemical potential difference is applied across the 
membrane. The bounds in these cases grow as iV 2 and N 1+l/ , respectively, where v is the exponent 
that relates the scaling of the radius of gyration to N. Our simulations demonstrate that the actual 
translocation times scale in the same manner as the bounds, although influenced by strong finite 
size effects which persist even for the longest polymers that we considered (N — 512). 

PACS numbers: 05.40.-a 05.40.Fb 36.20.Ey 87.15.Aa 



I. INTRODUCTION 



Translocation of a polymer through a narrow pore in 
a membrane is important to many biological processes, 
such as the injection of viral DNA into a host, DNA 
packing into a shell during viral replication, and gene 
swapping through bacterial pili [lj. Translocation also 
has practical applications in genetics as in cell transfor- 
mation by DNA electroporation , and in gene therapy 
Q ■ This has inspired a number of recent in vitro exper- 
iments, including the electric field-induced migration of 
DNA through microfabricated channels Q , or through an 
a-hemolysin protein channel in a membrane |fj . Ex- 
periments are motivated by the possibility to "read-off" 
a DNA or RNA sequence by tracking its passage through 
a pore 



Translocation of a polymer involves both molecular 
considerations, such as the shape of the pore channel 
and its interactions with DNA, as well as more macro- 
scopic factors such as the statistics and dynamics of the 
long polymer. It is the universal features of the latter 
which are the focus of this paper. While worming its way 
through the hole, the segments of the polymer on each 
side of the membrane can "explore" many possible config- 
urations. The number of allowed configurations actually 
is least when the polymer is halfway through the hole, 
presenting an entropic barrier. In this regard transloca- 
tion resembles other entropically controlled polymer sys- 
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FIG. 1: Schematic depiction of polymer translocation from 
the left side of a membrane to its right side. (Throughout 
the paper we shall follow this convention for the direction 
of transport.) The number s of the monomer at the hole is 
denoted the "translocation coordinate." 



terns, such as polymer trapping in random environments 
HUH, DNA gel electrophoresis [l^l or reptation [TT| . 
where the geometry of the obstacles constrains the kinet- 
ics of the polymer. 

A number of recent theoretical works have shed light 
on the translocation process jl 111 111 HI 111 111 fl7j . 
which is schematically depicted in Fig. ^ A single vari- 
able s representing the monomer number at the pore 
LLJ Ull HI] indicates how far the polymer has progressed, 
and is a natural variable for describing this problem. Due 
to its resemblance to the 'reaction coordinate' for chem- 
ical processes, we call s the translocation coordinate. If 
s changes very slowly, such that the polymer segments 
on both sides of the membrane have time to equilibrate, 
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the mean force acting on the monomer in the hole can 
be determined from a simple calculation of free energy, 
and the translocation problem is then reduced to the es- 
cape of a 'particle' (the translocation coordinate) over a 
potential barrier. In the following, we shall refer to this 
limit as Brownian translocation, but shall demonstrate 
that the required equilibration is not tenable for long 
enough polymers. 

In many experimental situations the polymers are not 
very long and the observed behavior strongly depends on 
the properties of the polymer and the pore. Thus intro- 
duction of some specific features into the channel proper- 
ties 0, 0] , or into interparticle potential [2(], EJ , may 
provide qualitative explanations of the observed behav- 
ior. While the theoretical understanding of the experi- 
ments is growing, we are still quite far from quantitative 
understanding of many observed features |22|. In this 
work, we restrict out attention to qualitative features of 
very long polymers that are independent of the details of 
the pore or inter-monomer potentials. Consequently, we 
restrict ourselves to simple models. Furthermore, since 
we are interested in contrasting the behavior of phantom 
polymers, in which the monomers do not interact with 
each other, with self-avoiding ones, in which monomers 
repel each other at short distances, we perform simula- 
tions in two-dimensional (2D) space, where such differ- 
ences are more pronounced. 

The reduction of the motion of a large molecule to 
a single-particle problem ignores the fact that the posi- 
tions of the monomers have strong correlations ^lj i lead- 
ing to non-trivial dynamical effects [Til |2^ | . In particu- 
lar, in a dilute solution of polymers in a good solvent, 
on time scales shorter than the overall relaxation time 
of the polymer, the motion of monomers is character- 
ized by anomalous dynamics p4l |25| . Not surprisingly, 
such effects are also present in translocation. In a pre- 
vious work we demonstrated [2(| that the scaling of the 
translocation times with the number of monomers has 
a power-law dependence which cannot be derived from 
Brownian motion of a particle over a barrier, but rather 
follows from general scaling considerations. In this work, 
we focus on consequences of anomalous dynamics in the 
presence of a force. 

The rest of this paper is organized as follows. In Scc.HTI 
we review the limit of very slow translocation, in which 
case the problem can be reduced to the Brownian motion 
of a single coordinate. We demonstrate that the require- 
ment of equilibration breaks down for long polymers, es- 
pecially in the presence of a force pulling the polymer to 
one side. We also demonstrate the importance of how 
such a force is applied to the polymer, contrasting the 
cases of a polymer pulled by the end, with one forced into 
a favorable environment. Lower bounds on the translo- 
cation times are obtained in Sec. lIIII bv comparison with 
the unimpeded motion of a polymer. The long time scale 
for equilibration of a forced polymer is due to its change 
of shape, e.g. into a stretched sequence of blobs when 
pulled at one end. Exactly how this change of shape 



is achieved through transmission of the force from one 
end to another is explicitly shown in Appendix IaI which 
solves this problem for a one-dimensional (ID) phantom 
polymer. Our other model system, the self-avoiding two 
dimensional polymer is discussed in Appendix [5] The 
bounds from unimpeded motion serve a convenient ref- 
erence point for discussion of anomalous processes that 
are described in detail in Sec. IIVI In particular, we find 
that the actual translocation times have the same scaling 
behavior as the lower bounds. 



II. BROWNIAN TRANSLOCATION AND ITS 
LIMITATIONS 

If the translocation process is sufficiently slow, at each 
stage the statistics of the segments will be governed by 
the equilibrium Boltzmann weight. If so, the dynam- 
ics is constrained to reproduce the corresponding statis- 
tics. A simple method to achieve this is to focus on the 
translocation coordinate s, and to write down a stochas- 
tic Langevin equation for its evolution. This procedure, 
which we shall refer to as Brownian translocation, is the 
chief tool employed in most analytical studies [lllll&ll8| . 
and shall be reviewed in this Section. 



A. Unforced motion 

In the absence of an external potential, the entropic 
contributions to the free energy of the two polymeric seg- 
ments result in a free energy J-'(s) — 7/csTln[(./V — s)s] 
[l2l [Til ITsI ] . Note that there is a decrease in the num- 
ber of possible states as the polymer threads the hole, 
which can be regarded as an entropic barrier. For phan- 
tom polymers (random walks) 7 is equal to 1/2, while 
for self-avoiding polymers it depends on the dimension- 
ality of space. From this free energy, we can construct 
a Langevin equation s = —mdTjds + rj(t), where m is 
a mobility coefficient indicating how easily the polymer 
is pulled through the hole. To obtain the correct Boltz- 
mann statistics, the noise ri(t) has to be uncorrelated at 
different times, with a variance equal to mksT. 

As demonstrated in Ref. [2(|, the corresponding 
Fokker-Planck equation for the probability p(s, t) can 
be made independent of N by a simple change of vari- 
ables s — ► sN and t — > t/N 2 . Consequently, the average 
translocation time r (and its fluctuations) must scale as 
N 2 for any 7. In fact, a very similar distribution for the 
transit time is obtained by ignoring the potential barrier 
completely (setting 7 = 0). In this limit the transloca- 
tion coordinate simply undergoes diffusion, i.e. at time 
scales much shorter than typical translocation times we 

expect {[s(t') - s{t)) 2 ^ oc \t - t'\. 
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B. Pulling on the ends 

The advent of optical tweezers has made it possible to 
manipulate single macromolecules. A common procedure 
is to attach latex balls to the end(s) of a polymer (such as 
DNAL and then to manipulate the balls by optical tweez- 
ers [27, 28] . While this is not the method commonly used 
in the translocation experiments, it motivates an interest- 
ing extension of the previous calculation. Let us imagine 
that through an optical tweezer set-up forces F + and F_ 
are applied to the two ends of the polymer (perpendicular 
to the wall), as in Fig. ^ 

A polymer configuration in which the ends are sep- 
arated by a distance r gets an additional Boltzmann 
weight of exp (^F ■ r/^rj . For a Gaussian polymer of 

length N, integration over all locations r leads to a con- 
tribution of exp [N(Fa/k B T) 2 ~\ to the partition function 
(a is a microscopic length scale). Restricting the inte- 
gration over r only to half space (as appropriate for the 
translocation problem), does not change the qualitative 
form, i.e., in the presence of the force the free energy 
acquires a term proportional to k B J 'N (F 'a / 'k B T) 2 . 

The above argument can be generalized to a self- 
avoiding polymer by noting that due to dimensional con- 
siderations a force F should always appear in the combi- 
nation FR/k B T, where R is a characteristic size of the 
polymer, such as its radius of gyration R g . Quite gen- 
erally R ~ aN v , with v = 1/2 for a phantom (non- 
self-interacting) polymer, and v — 3/4 and 0.59 for self- 
avoiding polymers in space dimensions d = 2 and 3, 
respectively [23. The contribution of the force to the 
free energy can thus be written as fegT$ {FaN u /k B T). 
When the force is sufficiently strong to deform the poly- 
mer into a sequence of blobs (we shall argue later that this 
is the relevant regime for translocation), the free energy 
is expected to be linear in N, necessitating ~ x 1 ^" ' . 
Such considerations thus imply a free energy contribution 
of Nk B T {Fa/k B T) 1,v . In the specific case of transloca- 
tion, adding up the contributions from the two segments 
leads to 



T{s) ~ k B T 



F + a 
kBT 



l/u 



(N- 



k^f 



l/u 



The corresponding Langevin equation for the translo- 
cation coordinate is now 



s = A 



(2) 



where we have absorbed various coefficients into the pa- 
rameter A. Note that the average velocity has a non- 
linear dependence on the forces; in the case of F_ = 
growing as FV V . Consequently, the translocation time 
in such a set up should decrease with the applied force as 
t oc N/F^ u . Note that this expression breaks down 
for small forces, where the typical translocation time 
is controlled by the diffusive fluctuations. The distinc- 



tion between weak and strong force regimes has a spe- 
cific meaning in the context of polymers, and quantified 
through the scaling combination / = FaN v /k B T . For 
weak forces this combination is small, and the equilib- 
rium polymer shape is not changed. For strong forces 
/ 3> 1, and the polymer becomes stretched. The same 
division applies to forces that are strong enough to over- 
come the diffusive character of the translocation coordi- 
nate. 



C. Chemical potential difference 

A more common situation for translocation is when 
the environments separated by the membrane are not 
equivalent, so that the polymer encounters a chemical 
potential difference between the two sides. In this case, 
the leading contribution to the free energy is T(s) = 
fj, + s + (N — s)/i_, and the Langevin equation takes the 
form 



s = A'A/j + 77(f), 



(3) 



with A/i = (/i_ — //_)_) > 0. In this case the average 
translocation velocity is predicted to be proportional to 
A/it, leading to typical exit times that scale as N/A/j, for 
large A/i. It is possible to envision situations in which the 
polymer is forced (or hindered) by a combination of both 
a chemical potential difference, and forces applied to the 
two ends. For the corresponding Langevin equation, we 
merely need to add the force contributions in Eqs. (0 
and 

In the experiments of Meller et al. |5j, translocation 
is driven by an electric potential difference between the 
two sides of an artificial membrane suspended in a liq- 
uid. It is commonly assumed that since the conductivity 
of the liquid is significantly higher than the conductivity 
of the membrane, the liquid on each side of the mem- 
brane is an equipotential. The voltage drop then occurs 
only across the membrane, and is experienced only by 
the (charged) monomers moving through the pore |22j . 
If so, the voltage difference is equivalent to the chemical 
potential difference discussed above. However, the previ- 
ous results with force acting only on the end monomers 
serve as a warning that the results are quite sensitive to 
where the force is applied to the polymer. It would thus 
be reassuring to carry out a more precise calculation of 
the electric field in the vicinity of the pore, and how it 
acts upon the monomers. 



D. Limitations 

The analytical procedure outlined in this section rests 
on the assumption that the two polymeric segments have 
come to equilibrium, so that the corresponding free en- 
ergy can be used to construct a Langevin equation. The 
minimal requirement is that the typical translocation 
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time r should exceed the equilibration time r cqu ii of a 
polymer. For a chain of finite size, it is always possible 
to achieve this by designing the pore to have a large fric- 
tion coefficient. In Ref. 15], it is argued that this is the 
case applicable to the experiments of Ref. |5j. However, 
the equilibration time of a polymer depends strongly on 
its length, scaling as r cqu ii(A) oc N zu . As discussed in 
the next section, the exponent zv is typically larger than 
2 for Brownian dynamics of self-avoiding polymers (and 
equal to 2 for a phantom polymer). 

Since typical transit times for unforced translocation 
scale as TV 2 , it is possible to imagine that the formalism 
may apply to phantom polymers. Indeed there is some 
evidence of this from numerical simulations 1^3] , although 
with an inexplicably large friction coefficient. The situ- 
ation becomes worse in the presence of a force (either 
imposed by pulling or a chemical potential difference) , in 
which case typical translocation times are predicted to 
be proportional to N. In the latter case, the range of 
applicability of 'Brownian translocation' is even further 
limited. 



III. UNIMPEDED MOTION OF A POLYMER 

Since the collective dynamics of the passage of polymer 
through the pore is hard to treat analytically, as a first 
step we shall derive lower bounds on the characteristic 
time scale. The key observation is that the restriction 
that the monomers should sequentially pass through a 
hole in a membrane can only impede the motion of the 
polymer. Hence the time scale for the polymer to travel 
the same distance in the absence of the wall should be a 
generous lower bound to its translocation time. In this 
Section we shall thus explore the time scale for unim- 
peded motion of the polymer in the circumstances of in- 
terest. 



A. Unforced diffusion 

In the unforced limit, the translocating polymer simply 
goes from one side of the membrane to the other by 'dif- 
fusion.' In the process, the center of mass of the polymer 
moves a distance of the order of a typical size, say the 
gyration radius of R g . How long does it take for a poly- 
mer to move a similar distance without the constraints 
imposed by the pore and the wall? In the absence of 
hydrodynamic interactions, the diffusion constant D of 
a polymer is related to the diffusion constant Dq of a 
single monomer by D = Dq/N. Consequently, the time 
that a polymer needs to diffuse its own radius of gyra- 
tion scales as N 1+2u [Tl| . (This is also the relaxation 
time of the slowest internal mode of a polymer , and 
is called the Rouse equilibration time.) For self-avoiding 
polymers v > 1/2, and the equilibration time is clearly 
longer than that obtained for translocation though a hole 
in the wall assuming Brownian translocation. The Rouse 



time scale of N 1+2v should thus be a lower bound to the 
correct translocation time. 



B. Pulling on the end 

Now consider a polymer that is being pulled by one 
end. The regime of interest to us is when the force is 
strong enough to deform the shape of the polymer. The 
equilibrium shape of the polymer is then a stretched se- 
quence of 'blobs' The number of monomers per blob 
Nb is such that force acting on it is marginally strong, 
and obtained from FaNg ~ fcgT. The typical size of 
each blob is thus Rb ~ aNg ~ ksT/F, while the number 
of blobs is N/N B ~ N(Fa/k B T) 1 /' J . The overall length 
of the stretched chain is now R(F) ~ Rb(N/Nb) ~ 
aNiFa/kBT) 1 /"- 1 . The mobility of the center of the 
mass of a polymer of length N is proportional to 1/N, 
and since there is only a force applied to one monomer, 
its net velocity scales as F/N . The characteristic relax- 
ation time is the same as the time scale of the polymer 
moving a distance of order of its size, and hence behaves 
as 

r cquil (F) ~ » oc N 2 F-^I». (4) 

(The same conclusion is obtained if we start with a glob- 
ular polymer and then apply the force to one end, and 
wait until the other end feels the force.) 

Note that upon approaching the boundary between 
weak and strong forces at F oc N~ w , we regain the equili- 
bration time N 1+2l/ for unforced polymers. However, the 
result in Eq. © is only valid for N~ v < Fa/k B T < 1, 
since for stronger forces, as we shall see in the specific 
models described in the next Section, the velocity of the 
monomer (and hence of the entire chain) saturates. 

Unimpeded motion of the pulled polymer thus places 
a lower bound of TV 2 on the translocation time, far ex- 
ceeding the time scale (oc N) calculated in the previous 
section. It may not be readily apparent how the force 
applied to one end of the polymer is transmitted to the 
other end, and why the qualitative picture of blobs pre- 
sented above is valid. There is actually one limit in which 
the problem of pulling a polymer by the end can be solved 
analytically, and that is for a ID phantom polymer. In 
Appendix [S] this model and the corresponding analysis 
are presented in some detail. 



C. Mimicking a chemical potential difference 

It is difficult to come up with an unhindered situation 
that best resembles the case of a chemical potential differ- 
ence across a membrane. Absent the restrictions imposed 
by the membrane, there is now a force that is applied 
to a single monomer, at the spatial position where the 
membrane would reside. Unlike the previous case, the 
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monomer to which the force is applied now changes con- 
stantly. There is thus no incentive for a drastic change in 
the shape of polymer, and we assume that the scaling of 
the size remains the same, i.e. R ~ aN v , independent of 
Afi. At each moment there is a force of A/z/a applied to 
the entire polymer, as a consequence of which its center 
of mass should move with a velocity v oc A/i/N. We thus 
conclude that the time for such unhindered motion scales 



, A . R N 1+u 

r(A/i) oc — — . 

V Afi 



(5) 



Note that to recover the equilibration time of the un- 
forced polymer we have to set A// <x N~ u in the above 
equation. While this is the same scaling form as that 
of a force applied to the end, it is different from the 
weak/strong criterion that would have been deduced on 
the basis of energetics (AfiN ~ fcsT). This is a reflec- 
tion of the manner in which we introduced the unimpeded 
version (as a force, rather than a chemical potential dif- 
ference.) Nonetheless, we still expect the velocity, and 
hence the time scale in Eq. (JjjJ to saturate for Afi ~ fc^T, 
as explicitly demonstrated for the models considered in 
the next Section. 
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FIG. 2: Temporal fluctuations of the monomer number s 
located at the hole, averaged over the initial time t and over 
1000 independent simulations, for polymers of lengths N = 8, 
16, 32, 64, and 128. 



suggested 26] that the short time fluctuations of s(t) fol- 
low the anomalous diffusion relation 



IV. ANOMALOUS TRANSLOCATION 

Having established some some (presumably generous) 
lower bounds, we now would like to focus on the true 
asymptotic dynamics of translocation. Given the limita- 
tions of analytical studies, the chief tool employed in this 
section is numerical simulations. Interestingly, we find 
that the lower bounds obtained in the previous section 
are actually quite restrictive. 

A. Sub— diffusive behavior of unforced motion 

In a previous work p6| , we made a detailed study of the 
iV-dependence of the mean translocation time. One of 
the central conclusions was that in the case of Brownian 
dynamics of a self-avoiding polymer, the translocation 
time scales as 

r ~ N 1+2 ». (6) 

This is of the same order as the equilibration time of 
a polymer of length N, and also demonstrates that the 
actual exit time scales in the same manner as the bound 
established in the previous section. 

The above A-dependence of r is inconsistent with sim- 
ple diffusion of the translocation coordinate s(t), reflect- 
ing the constraints imposed by the collective motion of 
the entire polymer. A related situation occurs for the 
fluctuations of a labeled monomer in space, which are also 
anomalous and sub-diffusive p4j on time scales shorter 
than the equilibration time. Following this analogy, we 



{As 2 (t))~t^. (7) 

For Eqs. © and Q to be consistent, we must obtain 
As of order N, when t of order r, leading to the expo- 
nent relation ( = 1/(1 + 2v). Note that for a phantom 
polymer v = £ = 1/2, i.e. the anomaly disappears in 
this limit, and the process becomes diffusive. (This dif- 
fers from the corresponding motion of a labeled monomer 
|24l |25| , which remains anomalous even for a phantom 
polymer.) This is consistent with a detailed study of a 
three-dimensional phantom polymer by Chern et al. |30l | 
which concluded that the results may be interpreted in 
terms of diffusive motion of the translocation coordinate 
over a barrier. Such correspondence is likely a fortuitous 
coincidence for phantom polymers, and even in this case, 
the value of the effective diffusion constant could not be 
obtained from the geometrical features of the model [3(J • 
To establish the anomalous nature of (unforced) 
translocation dynamics, we carried out Monte Carlo 
(MC) simulations on a model of two-dimensional self- 
avoiding polymers, described in Appendix [5] Simula- 
tions in two (rather than three) dimensions have the ad- 
vantage of relative ease, and stronger differences from 
phantom polymers. We followed the dynamics of s(t), 
focusing on the quantity As 2 = (s(t' + 1) — s(i')) 2 . This 
correlation function is depicted in Fig. |2 and was ob- 
tained by averaging over t' of over 1000 independent sim- 
ulations for N — 8, 16, 32, 64 and 128. Two cautionary 
points must be made in considering this data: The first 
is that we have no a priori assurance that this process is 
stationary; the results may depend on both t and t', and 
consequently influenced by the averaging over t' . (For 
example, choosing a short averaging range may increase 
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FIG. 3: "Snapshots" of a 128-monomer phantom polymer 
passing through a membrane in one dimension. Each line 
depicts the position x (in lattice units) of the nth monomer 
at a fixed time. Different lines correspond to times when 
the 20th, 40th, . . . 120th (bottom-left to top-right) monomer 
crosses the membrane (the thick line at x — 0). 



FIG. 4: Logarithmic plot of the mean translocation time as 
a function of length for a one-dimensional phantom polymer 
with an infinite force applied to one end. The circles repre- 
sent passage through an opening, while the squares represent 
motion in the absence of a membrane. Each data point rep- 
resents an average over 1000 processes. 



the effect of the initial conditions. We performed aver- 
aging over shorter ranges of t' and did not see signifi- 
cant differences in the correlation functions.) Secondly, 
for each value of t, we can only include processes whose 
translocation time exceeds t\ consequently, the size of the 
ensemble decreases with increasing t. (However, this ef- 
fect is insignificant for t several times shorter than the 
mean translocation time.) In Fig. [3 all available t' are 
included, i.e. for each polymer the statistics was collected 
up to the moment that the translocation was completed. 

Since the values of s cannot exceed N, as the time- 
difference t becomes of order of t, As 2 saturates, as is 
apparent in the case of N = 8 in Fig. However, for 
times shorter than t(N), the results for different lengths 
seem to form a single curve. This instills confidence in the 
quasi-stationary character of translocation on time scales 
shorter than r. On the logarithmic scale the graph seems 
to have curvature for t < 1000. This is probably a conse- 
quence of discreteness of the model, since corresponding 
differences in s are smaller than 5. For larger times the 
slope of the curve approaches 0.80, which clearly indi- 
cates the presence of anomalous diffusion, and is consis- 
tent with the expected value 2/(1 + 2z/), with v = 3/4 
for 2D self-avoiding walks [29j. Thus, despite its finite 
duration, the translocation process at short scales resem- 
bles a stationary process (at least As 2 is insensitive to 
t') which exhibits anomalous dynamics. 



B. Pulling on the end 

As discussed in Sec. [H] the polymer pulled by a 
force F + is relatively undistorted as long as / = 



F + aN v /JcbT <C 1, and the corresponding translocation 
times are not very different to those in the absence of 
force. Increasing polymer length at fixed F + ultimately 
leads to a regime with / 3> 1, in which the polymer is 
expected to be stretched into a sequence of blobs. It is 
the latter regime which is of interest us, and which shall 
be explored by examining the one- and two-dimensional 
polymer models introduced earlier. 



1. One-dimensional phantom polymer 

Simulations are carried out with the model ID phan- 
tom polymer presented in Appendix ^ starting with a 
polymer that is equilibrated on one (say, left) side of 
the "membrane" (a point on a ID lattice) with one end 
point held at the membrane. The "narrow opening" is 
implemented by allowing only sequential passage of the 
monomers across the membrane, i.e. the nth monomer 
can move from the left to the right only if (n — l)th 
monomer is already on the right side. Conversely, the 
nth monomer may diffuse from the right side to the left, 
only if the (n + l)th monomer is on the left side. The 
first monomer of the chain is restricted to remain on the 
right throughout the process. We study the dynamics of 
translocation as a function of the force F applied to the 
first monomer. However, the results become independent 
of F when the reduced force / = Fa/kBT exceeds unity, 
since it becomes very unlikely for the first monomer to 
move backwards. 

Figure [31 depicts a sequence of "snapshots" of a 128- 
monomer polymer going across the membrane. Since the 
maximal separation between adjacent monomers is two 
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lattice spacings, the slopes of the curves are limited by 
2. Note that for x > the polymer is almost maximally 
stretched, while the x < configurations resemble the 
initial random walk state. There is also much similarity 
between the profile of the polymer for x > 0, and the 
the steady-state configurations of a polymer moving in 
the absence of the membrane, as depicted in Fig. [21 of 
Appendix [X] 

The results of averaging the translocation time (over 
1000 realizations) are indicated by the circles in Fig. 0] 
The points appear to fall on a straight line in this loga- 
rithmic plot, with the slope of 1.93 ± 0.01 from a least- 
squares-fit. There is a slight upwards curvature, and the 
effective slope varies from 1.84 for points with N < 128 
to 1.93 for all the points, indicating potential crossover 
effects persisting even for N = 512. From this data by 
itself it is difficult to determine the ultimate slope. How- 
ever, we can compare the results with the times required 
to cross an imaginary membrane (i.e. unimpeded diffu- 
sion). This lower bound which was described in the pre- 
vious section (and discussed in detail in Appendix |A"|) . 
leads to the mean passage times depicted by the squares 
in Fig.Q] The extrapolated slope for this unimpeded mo- 
tion is indeed 2.00. Since the unimpeded crossing times 
are indeed shorter, the asymptotic exponent in the pres- 
ence of the membrane has to be larger than or equal to 
2. (Otherwise, for sufficiently large N the curves will in- 
tersect causing longer times for passage if the membrane 
is absent.) We therefore conclude that the translocation 
of a phantom polymer in ID should asymptotically scale 
as N 2 , saturating the bound obtained previously. 



2. Two-dimensional self-avoiding polymer 

We next study the translocation of a self-avoiding 
polymer in 2D as a function of the force F applied to 
the first monomer. Figure [5] depicts the distribution of 
translocation times for a 128-monomer polymer at three 
values of / = Fa/ksT. As / increases from 0.25 to oo, 
the mean translocation time drops by less than one order 
of magnitude, and the relative width of the distribution 
decreases somewhat. Note that once / 3> 1, the first 
monomer always moves in the forward direction, and the 
results become independent of /. 

Figure [S] summarizes the results obtained for poly- 
mer lengths N ranging from 8 to 128, and for a vari- 
ety of forces. Each point corresponds to an average over 
1000 realizations. The figure depicts the scaled inverse 
translocation time as a function of the dimensionless force 
/ = Fa/ksT. The vertical scale has been multiplied by 
N 7 to produce moderate collapse of data for different 
Ns, although as explained further on, we do not believe 
this to be the correct scaling factor in the asymptotic 
regime. As expected, the curves saturate when / signifi- 
cantly exceeds unity. 

Note that all the points in Fig. belong to the regime 
where / 3> 1, i.e. when the shape of the polymer is ex- 
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FIG. 5: The distribution of translocation times for a 2D 
polymer with N = 128. Each histogram represents results 
from 250 independent translocations for forces FajksT = 
0.25, 1, and oo (from right to left) applied to one end. (The 
horizontal axis is logarithmic.) 
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FIG. 6: Scaled inverse translocation time as a function of 
the reduced force Fa/ksT applied to the end-monomer, for 
AT = 8, 16, ...128. 



pected to be different from equilibrium, and stretched. 
This is confirmed in Fig. which depicts configurations 
of the polymer in the process of translocation under the 
action of an infinite force. The front end of the poly- 
mer is quite stretched, somewhat resembling a directed 
random walk, suggesting that self-avoiding interactions 
play a secondary role in this limit. If the front part of the 
polymer controls the translocation time, it should have 
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FIG. 7: Configurations of a polymer of length N = 128, 
pulled through a hole by an infinite force applied to its first 
monomer. The circles, diamonds and triangles represent the 
initial configuration, and at times t = 60, 000 and 120,000 
Monte Carlo time steps. 
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FIG. 8: Logarithmic plot of the dependence of the translo- 
cation time r on polymer length N, when an infinite force is 
applied to the end-monomer. The line is a fit to a power law 
dependence with exponent 1.875. 



the same scaling with N as the corresponding time for 
a phantom polymer, i.e. we expect r oc A^ 2 . Figure 
depicts the dependence of r on N; the effective exponent 
in this range is 1.875 ± 0.005. Although smaller than 
two, it is close to the value of the effective exponent of 
a phantom polymer in this range of lengths N. Given 
the bound presented earlier, it is reasonable to expect 
asymptotic convergence to this value. However, as far as 
we can judge by analogy to phantom polymers, we need 
N to be much larger than 1000 to see an exponent of 2. 



C. Chemical potential difference 

When the environments on the two sides of the mem- 
brane are different, the monomer at the pore experiences 
a force pushing it to the more favorable side. As ex- 
plained previously, this form of forcing leads to yet a 
different form of asymptotic behavior which is once more 
explored using our two numerical models. 
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FIG. 9: Logarithmic plot of the mean translocation time as 
a function of polymer length for a one-dimensional phantom 
polymer subject to an infinite chemical potential difference. 
The solid line is a fit to a power law with exponent 1.45. 

1. One-dimensional phantom polymer 

Figure El depicts the A-dependence of the transloca- 
tion time t for a ID phantom polymer under the influence 
of an infinite potential difference, i.e. when the monomer 
at the pore can only move to one side. The data on the 
logarithmic plot are fitted to a straight line with exponent 
1.45 ± 0.01, although there is a slight upward curvature 
even for A~ = 512. Note that in the limit of large chemi- 
cal potential difference, the monomers that have already 
crossed to the right side no longer play any role in the 
translocation process, which is thus constrained by the 
dynamics of the monomers remaining on the left side. 
In Sec. 11111 we argued that the translocation time of an 
unimpeded phantom polymer should scale as A~ 3 / 2 , since 
the leftmost monomer must travel a distance of order 
N 1 / 2 with a velocity of order 1/JV. Thus 3/2 should be 
a lower bound for the exponent characterizing the scal- 
ing of the translocation time with N. By comparing this 
limit with our numerical results, we conclude 3/2 to be 
the true asymptotic form describing our simulations. 

2. Two-dimensional self-avoiding polymer 

Figure [23 presents the distribution of translocation 
times for a polymer with A^ = 64 at several values of 
the dimensionless chemical potential difference A^/fceT. 
These distributions are quite wide — although they be- 
come (relatively) narrower with increasing A/j,, the width 
of the distribution is of the same order as the average even 
in the limit of an infinite A/x. As expected, the average 
translocation times decrease and saturate when Afi/kgT 
exceeds unity. The results for different values of A/i/hsT 
and A~ can be approximately collapsed by scaling r with 
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FIG. 10: Distribution of translocation times of a 64- 
monomer polymer subject to chemical potential differences 
of An/k B T = 0, 0.25, 0.75, and 2 (right to left). (The hori- 
zontal axis is logarithmic.) 
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FIG. 11: Scaled inverse mean translocation time as a function 
of reduced chemical potential difference A/i/fcsT, for N — 8, 
16, 32, 64 and 128. 



TV 1 - 45 , as shown in Fig. 1111 (where each point represents 
an average over 1000 independent runs). The quality of 
the collapse is very poor, and we shall argue that iV 1 ' 45 
is not the expected asymptotic power. 

Since (with the exception of A/i = 0), the points in 
Fig. II II correspond to a strong force at the pore, we ex- 
pect the configurations of the translocating polymer to 
be different from those of a polymer in equilibrium. To 
explore this difference, in Fig. I12l we show a pair of con- 
figurations for an infinite chemical potential difference. 
We see that fast translocation results in a higher den- 
sity of monomers immediately to the right of the pore, 
which may in principal slow down the process. (Recall 
that in the case of phantom polymers the monomers that 




FIG. 12: Configurations of a polymer of length N = 64 cross- 
ing a membrane from left to right under an infinite chemical 
potential difference. Full and open symbols represent times 
of t = 10, 000 and 25,000 Monte Carlo steps, respectively 



have passed through the hole have no further influence.) 
Nevertheless, the whole process should still be bounded 
by the corresponding time for passage of an unimpeded 
polymer, as discussed in Sec. IIIII For considerations of 
this bound, the relevant time corresponded to motion of 
the leftmost monomer over a distance of size R oc N", 
leading to a time scale growing as N to a power of 1 + v 
which in 2D is 1.75. 



The optimal data collapse of the inverse translocation 
times for N < 128 leads to an exponent 1.45 (with a 
rather poor data collapse). Therefore, we extended our 
simulations for the case of the infinite A/i to larger values 
of N. Fig.^|]represents the dependence of mean translo- 
cation time on TV, for polymer lengths up to 512. Data 
points with N < 128 correspond to averages over 1000 
independent simulations, while N — 256 and 512 include 
300 and 130 runs, respectively. The effective slope of the 
fit for data points below N = 128 is 1.45 ±0.01, while all 
the data points produce an effective slope of 1.53 ± 0.01. 
Also by directly measuring the effective slopes between 
successive pairs of points we definitely see an increase, 
with the last pair of points giving a slope of 1.60 ± 0.03. 
However, the increase is very slow, and the uncertainties 
are too large to enable a reliable extrapolation to large 
N. By comparing the results to data obtained in the sim- 
ulations of phantom polymers we believe that eventually 
the exponent will reach 1.75; however, this will probably 
happen only for N significantly larger that 1000. 
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FIG. 13: Logarithmic plot of the dependence of the translo- 
cation time t on polymer length N, for an infinite chemical 
potential difference. The solid line is a fit to a power law 
dependence with exponent 1.53. 



V. DISCUSSION 

Translocation of a polymer through a pore is intrinsi- 
cally a many body problem involving collective and coop- 
erative motion of monomers crossing a membrane. If the 
process is sufficiently slow, it is possible for the segments 
on the two sides to come to equilibrium, in which case the 
dynamics of the translocation coordinate (the number of 
monomers on one side) is similar to Brownian motion of 
a single particle. However, the assumption of equilib- 
rium must necessarily break down for long enough poly- 
mers; even more drastically when the polymer is pulled 
to one side by a force. We argued previously [2|| that 
the collective motion of the whole polymer slows down 
the translocation to the extent that the corresponding 
dynamics is anomalous and sub-diffusive. This is explic- 
itly verified in this paper by quantifying the temporal 
correlations of the fluctuating translocation coordinate. 
We note that there is a general theoretical framework |3l| 
for anomalous dynamics of a single variable, which may 
be profitably applied to this problem. 

One of the objectives of this paper was to find how dy- 
namical anomalies affect the motion of the polymer under 
the action of a force. With this in mind, we also empha- 
sized that the method by which the polymer is forced 
is quite important. In particular, pulling the polymer 
by one end leads to stretched configurations, and slower 
overall dynamics, compared to applying a chemical po- 
tential difference (which can modify the densities on the 
two sides). While pulling the polymer by optical tweezers 
is not the currently favored method for artificial translo- 
cation of biopolymers, for potential applications such as 



decoding the sequence, it should offer a better controlled 
procedure (whether by itself or in conjunction with a 
voltage difference). 

To understand the time scales involved in forced 
translocation, we initially provided what at first glance 
appear to be a quite loose lower bounds by analogy to 
unimpeded motion of a polymer, i.e. neglecting the con- 
straints imposed by passage of the polymer through a 
hole in a wall. We then performed numerical simulations 
on two model systems: a ID phantom polymer, and a 
self-avoiding polymer in 2D. Direct interpretation of the 
numerical results was made difficult by very large cross- 
over effects which persist in the length scales of 100-1000 
monomers accessible to numerical study. Nonetheless, 
by appealing to the lower bounds found earlier, we con- 
cluded that (rather surprisingly) the actual translocation 
times scale in the same way as in the limit of unimpeded 
motion. Thus the constraints from the collective motion 
of the whole polymer turn out to be at least as important 
as those imposed by the requirement of passage through 
a hole. 

The experiments of Ref. [j| suggest that in the range 
of 10 to 100 base pairs, the pulling velocity of single- 
stranded DNA through a nanopore is independent of N, 
but with a non-linear dependence on the applied force. 
In Ref. |2y] , we briefly speculated whether such behavior 
may be consistent with anomalies associated with poly- 
meric constraints. As demonstrated in this paper, such 
constraints result in time scales (and hence pulling veloc- 
ities) which must depend on N for large enough N. The 
only case where we observe a non-linear force-velocity 
relation which is independent of N is when short (hence 
equilibrated) polymers are pulled by a force applied to 
one end. 

We hope that these results encourage further experi- 
mental and analytical studies of forced translocation. In 
particular, it would be interesting to better characterize 
the manner in which external forces act on the polymer, 
even in the case of a voltage difference across the mem- 
brane. Hydrodynamic effects, not considered in this pa- 
per, are also likely to play an important role. Some of 
these effects can be included in more realistic numeri- 
cal simulations, although in that case one should keep in 
mind the rather long crossover times that appear to be 
intrinsic to this process. 
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APPENDIX A: THE ONE-DIMENSIONAL 
PHANTOM POLYMER MODEL 

Self-avoiding interactions are integral to understand- 
ing the statics and dynamics of real polymers. Never- 
theless, it is useful to study one-dimensional phantom 
polymers, with no interactions between monomers which 
are not adjacent along the chain. Independently of the 
details of the interactions of adjacent monomers, long 
phantom polymers are harmonic |Tlj . in the sense that 
the probability distribution of the distance |r$ — | be- 
tween monomers i and j, approaches a Gaussian for large 
\i — j\. In practice, such behavior already appears for 

— ~ 10, and on such a "coarse-grained" level one can 
view the polymer as consisting of monomers connected 
by springs whose energy is proportional to (r^ — rj) . In 
this limit, certain aspects of phantom polymer dynamics 
can be analyzed analytically, and we can compare the 
expected asymptotic behavior with the numerically ob- 
served dynamics. Such treatment provides both a better 
view of crossover effects, and produces some insights into 
the dynamics of more realistic (self-avoiding) models. 

We employ a phantom polymer model in which the 
monomers are restricted to sites of a ID lattice with spac- 
ing a. The polymer connectivity is implemented by re- 
quiring the distance between adjacent monomers not to 
exceed two lattice constants. An elementary Monte Carlo 
(MC) step consists of randomly picking a monomer and 
attempting to move it in a randomly selected direction. If 
an external force F is applied to the first monomer of the 
chain, then the probability to make a step in the direction 
opposite to the force is proportional to exp(—aF/ fc.eT) . 
A MC time unit corresponds to N attempts to move 
monomers. 

We first performed simulations of polymer motion 
when F — > oo, in the absence of a membrane. Fig- 
ure El depicts several examples of spatial configurations 
in steady state, i.e. when the initial spatial configuration 
has been forgotten. The resulting profiles can be ex- 
plained analytically by noting that in the Gaussian limit 
the equation of motion of a ID phantom polymer is given 
by 

< Ai > 

where we have used dimensionless units, (distance in lat- 
tice constants, and time in MC units, which leads to a 
(monomer) diffusion constant of order one) , and omitted 
prefactors of order unity. The model used in the MC 
simulations has a finite stretchability, while a Gaussian 
chain can be stretched indefinitely. However, once the 
scaled force / = Fa/ksT becomes significantly larger 
than unity, the dynamics becomes essentially indepen- 
dent of F, as a backwards step of the end point has 
negligible probability. Thus, instead of using F — » oo, 
we can get set / = 1, which leads to the boundary con- 
dition dx/dn = — 1 at the beginning of the chain. For 
such a boundary condition we can easily find a stationary 
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FIG. 14: Five profiles of a ID phantom polymer of length 
N = 128, with an infinite force applied to one end in steady 
state motion. The curves are displaced along the vertical axis 
such that the position of the first monomer is at x = 0. 



2 




1 2 3 4 5 

t/N 2 

FIG. 15: Scaled coordinate x of the first monomer as a func- 
tion of scaled time t, for N = 8, 16, ... , 512 (right to left), in 
the absence of a membrane. 

solution of Eq. i|Al|) , as 

x{n) = vt+^(n~N) 2 , (A2) 

where v = 1/N. Simulations, of our model polymer in 
the absence of a membrane indeed confirm that the sta- 
tionary state velocity is proportional to 1/N. Moreover, 
the actual shape of the polymer in steady state, as de- 
picted in Fig. 1141 is similar on average to Eq. JHJ|. 

The solution in Eq. (|A1(1 can also be used to under- 
stand the crossover to stationary motion, starting from 
a relaxed initial state. It is easy to see that the initial 
velocity of a point to which the force is applied is of order 
unity. As time progresses, this velocity decays as i -1 / 2 , 
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of the motion of the first coordinates should collapse for 
different values of N. The actual picture appearing in 
Fig. El is more complicated: while we see an approach 
to a single curve for the largest N, there are very strong 
finite size effects for moderate values of N. The results 
for crossing the imaginary wall are surprisingly similar to 
those for translocation of the phantom polymer through 
a hole (as described in Sec. IIVf) . as depicted in Fig. I16I 
(Both plots are obtained by averaging 1000 translocation 
processes.) The linear behavior of these curves close to 
the origin confirms our expectation that at short times 
the velocity is proportional to t~ x l 2 . However, for scaled 
variables around 0.5 the line has a slight curvature, which 
distorts the apparent scaling relations, and creates the 
illusion of slightly different exponents. 



FIG. 16: Scaled coordinate x of the first monomer of a 
translocating phantom ID polymer, as a function of the scaled 
time t, for N = 8, 16,. . . , 512 (right to left). 

until after time r w N 2 it reaches its final (stationary) 
value of order 1/N. The monomer at the opposite end 
does not feel the external force in the beginning, and 
starts moving with velocity of order 1/N after time t. 
We also performed simulations of unimpeded polymer 
motion mimicking the translocation set-up, by consider- 
ing an imaginary membrane located at position x = 
which has no effect on the motion of monomers. The ini- 
tial configuration was chosen by equilibrating the poly- 
mer on one side (x < 0), with the first monomer fixed 
to x = 0. Then an infinite force was applied to the first 
monomer, and its position x was tracked as a function of 
time, until all the monomers crossed to x > 0. For every 
N, the results were averaged over 1000 independent runs. 
In these simulations the end point only needs to move a 
distance of order R g ~ N 1 / 2 to cross to x > 0. Moving at 
a steady velocity of 1/N, this would take a time of order 
iV 3 / 2 , which is significantly shorter than r. Thus the time 
~ -ZV 2 required for the last monomer to start feeling the 
force sets the time scale for a polymer with a large force 
applied to its first monomer, to move a distance of order 
of its radius of gyration. The measured iV-dependence 
of the translocation time is depicted by squares in Fig. ^ 
There is a slight curvature in the logarithmic plot and 
the slope approaches 2.00 ± 0.01, confirming r ~ TV 2 . 

The above arguments suggest that (x/N) 2 should be a 
linear function of t/N 2 . Thus appropriately scaled plots 



APPENDIX B: THE TWO-DIMENSIONAL 
SELF-AVOIDING POLYMER MODEL 

Wc used a 2D lattice fluctuating bond polymer model 
p^j for MC simulations of a self-avoiding polymer. The 
monomers are placed on the sites of a square lattice, with 
the bonds between adjacent monomers restricted not to 
exceed ylQ lattice constants. The excluded volume be- 
tween monomers is implemented by requiring that no 
two monomers can approach closer than 2 lattice con- 
stants. The membrane with a hole is constructed from 
a row of immobile monomers arranged in a straight line, 
with a 3 lattice constant gap representing the hole. Such 
a hole is small enough to allow only a single monomer 
to pass through, thus enabling a unique identification of 
the monomer s which separates the polymer into two seg- 
ments on different sides of the membrane. An elementary 
MC move consists of randomly selecting a monomer and 
attempting to move it onto an adjacent lattice site (in 
a randomly selected direction). If the new position does 
not violate the excluded-volume or maximal bond-length 
restrictions, the move is performed. N elementary moves 
form one MC time unit. The first monomer is not allowed 
to withdraw to the opposite side of the membrane. Since 
we are investigating a non-equilibrium process, the ini- 
tial conditions may play an important role. We chose an 
initial state in which the first monomer was fixed inside 
the hole, and the remaining polymer was equilibrated for 
more than the Rouse relaxation time After such 

equilibration, the first monomer was released, and that 
moment was designated as t = 0. 
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